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Abstract 

In the (2+l)-dimensional QED with and without the Chern-Simons term, we find 
the non-local gauge in which there is no wavefunction renormalization for the fermion 
in the framework of the Schwinger-Dyson equation. By solving the Schwinger-Dyson 
equation in the non-local gauge, we find a finite critical value Nj for the number of 
flavors Nf of four-component Dirac fermions, above which the chiral symmetry restores 
irrespective of presence or absence of the Chern-Simons term. In the same framework, 
we study the possibility of dynamical breakdown of the parity. It is shown that the 
parity is not dynamically broken. We discuss this reason from the viewpoint of the 
Schwinger-Dyson equation. 
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1 Introduction 



The study of the dynamical symmetry breaking in the (2+l)-dimensional QED (QED3) 
has attracted a great deal of attention recently. In the framework of the Schwinger- 
Dyson (SD) equation in the relativistic quantum field theory, Appelquist et al. fl], 
investigated the dynamical symmetry breaking of U(2Nf) to U(Nf) x U(Nf) in QED3 
with Nf flavors of four-component Dirac fermions. The study has given rise to much 
controversy as follows. In the bare vertex approximation: T^(p,q) = 7„, Appelquist, 
Nash and Wijewardhana (ANW) have shown that there exists the finite critical 
number of flavors 

Nf=-z~ 3.2, (1) 

above which (Nf > Nf) the chiral symmetry restores Q and for Nf < Nf the fermion 
mass is dynamically generated and obeys the scaling of the essential singularity type: 



m ~ Af(Nf), f(N f ) = exp 



2nir 



^jNf/Nf - 1 



(2) 



This result was obtained under the assumption of no wavefunction renormalization 
A(p) = 1 for the fermion in the Landau gauge where the fermion propagator is written 
as S(p) = [A(p)f— B(p)]~ 1 . They justified their approximations based on the argument 
of l/Nf expansion in which A(p) = 1 + 0(l/N } ) and T^p, g) = 7(U + 0(1/N f ). 

However, this was criticised by Pennington and Webb (PW) [|]] and Atkinson, 
Johnson and Pennington (AJP) ||. They claimed that the one loop correction to 
the wavefunction renormalization and the corresponding vertex correction of the type 
r^(p, q) = l^A(q) (] leads to the scaling law of the exponential type: 



f(N f ) = exp 



37T 2 



-N 



(3) 



and that there does not exist a finite critical flavour, Nf = oo. Similar claim was 
stated in the approach of the effective potential in contrast to the original result by 
Matsuki, Miao and Viswanathan [[?J. 

On the other hand, the Monte Carlo simulation of non-compact QED on a lattice by 
Daggoto, Kogut and Kocic Q suggests the existence of two phases which are separated 
by a finite NS ~ 3.5 ± 0.5. However Azcoiti et al. Q recently report different and 
delicate results. Hence this issue has not yet been confirmed by the simulation of 
lattice non-compact QED. 

In the framework of the SD equation, it was also pointed out that the inclusion of the 
infrared (IR) cutoff which plays the role of the finite size of the lattice simulation may 
drastically change the critical phenomenon associated with the dynamical symmetry 
breaking. Especially the ansatz of PW leads to the mean-field type scaling law and the 



1 In the quenched limit Nf — > 0, the theory has only one phase, chiral-symmetry-breaking phase, as shown 
by the SD equation and the Monte Carlo simulation ||] . 

2 Similar ansatz for the vertex, e.g., T^p, q) = ±7 M [A(p) + A(q)] or T^faq) = "f^[A(p)6(p 2 - q 2 ) + 
A(q)9(q 2 — p 2 )}, leads to the same result. 
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cutoff-dependent finite Nf in the large IR cutoff, while the exponential type scaling 
and infinite are obtained in the sufficiently small IR cutoff |pX| l. 

The subtlety of the problem in the framework of the SD equation comes from 
the fact that, in order to take into account the wavefunction renormalization, it is 
indispensable to include the vertex correction as a consequence of the Ward-Takahashi 
(WT) identity. However incorporating the vertex correction properly is quite difficult 
in the non-perturbative study of gauge theory, although there are many works of going 
along such a line under the appropriate ansatz for the vertex in QED3 [11, 12 1 as well 
as the (3+l)-dimensional case, QED4 

In this paper we investigate this problem from a different point of view. Instead of 
considering the vertex correction, we look for the situation in which we do not have to 
take into account the vertex correction. This is achieved only when there is no wave- 
function renormalization. In the quenched approximation which neglects the vacuum 
polarization in the photon propagator, this is simply realized by taking the Landau 
gauge @. In the presence of the vacuum polarization and/or the Chern-Simons (CS) 
term [|l7|], however, we cannot choose the (covariant) gauge fixing parameter such that 
A(p) = 1 follows. Recently it has been recognized that this can be done by taking the 
non-local gauge 



ia 19 



|20|| . This standpoint leads to the complicated non-local gauge 
instead of the vertex correction. However this method has some advantages that we 
do not have to treat the coupled SD equation for two functions A(p) and B(p), and 
somewhat arbitrary (non-perturbative) ansatz for the vertex so as to satisfy the WT 
identity. 

In the same scheme, we study the dynamical breakdown of the parity. The bare 
CS term breaks the parity explicitly. Even if there exists no bare CS term, the CS 
term may be induced by the radiative correction through the non-perturbative effect. 
Therefore there is a possibility that the induced CS term breaks the parity. However 
previous studies so far support the claim that the dynamical breakdown of the parity 
does not occur [21, 22, 23, in agreement with the general argument 2Bj. We clarify 
this reason from the viewpoint of SD equation. 



2 SD equation and the non-local gauge 

We consider the (2+l)-dimensional QED (QED3) with the following (euclidean) la- 
grangian 

£qed3 = <F(«7/A -m e - m T)i^ + ~F* V + eft-y^A^ + £ C s + £gf, (4) 
with the Chern-Simons term 

i 

£cs = ^e^ up A^d u A p , (5) 

and the gauge fixing term Cqf whose form will be specified in what follows. According 
to Appelquist et al. Q, ip l denotes the 4-component Dirac fermion with a flavour index 
i(i = 1, ...,Nf). Here 7 a (a = 0, 1, 2) are 4x4 matrices which satisfy the Clifford algebra 
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{linlv} = —l&iiv an d are explicitly defined by [] 

f-ia 3 \ fia 1 \ fia 2 \ . 

with o" a (a = 1, 2, 3) being the Pauli matrices. We further introduce the 4x4 matrices 
73, 75 and r as 

73:= (i o) ' 75 := 70717273 = (i o 1 ) ' r := 7375 = (o -l)' ^ 

First we define the chiral transformation by ip(x) — > 73^(2;), A^(x) — > A^x). 
Then the mass term me^tp and m ^Tijj are respectively odd and even under the 
chiral transformation. Next the parity transformation is defined by ^(xq, X\, x 2 ) — ¥ 
Pip(xo,-xi,X2) 2 A fl (xo,x 1 ,x 2 ) -» (-lY^A^xo, -xi,x 2 ), with P = -27571. Then 
m e i^ijj and iJioipTip are respectively parity even and odd mass term. The bare CS 
term Cos is parity odd and hence the inclusion of the bare CS term breaks the parity 
explicitly even if m = 0. 

Corresponding to the diagram of Fig.l, the SD equation for the full (exact) fermion 
propagator S(p) = [A{p)f — B{p)]^ 1 is written as follows. 

A{p)i> — B(p) = p 1 — m e — m T 

1 2 f d 3 q A{q)i-B{q) 

where D fJiU (k) is the full photon propagator and T u (p,q) the full vertex function. In 
this paper we take the bare vertex 7^ instead of the full vertex r«(p, q). In order for the 
bare vertex approximation T^p, q) = 7^ to be justified in this SD equation, we must 
guarantee no wavefunction renormalization A(p) = 1 in light of the Ward-Takahashi 
(WT) identity: (p„ - q^T^q) = S" 1 ^) - S^fa). Actually, in the case of the 
quenched approximation (without CS term): 
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(9) 



for the covariant gauge-fixing Cgf = ^(c^A^) 2 , this is simply achieved by taking the 
Landau gauge £ = 0, see e.g. ||. 

In the presence of the vacuum polarization and/or the CS term, this cannot be 
achieved by choosing an appropriate value for the usual gauge fixing parameter £ and 
instead we must adopt the non-local gauge in the sense that £ becomes a function of 
the photon momentum: £ = £(/c 2 ) ||, 0, [Hj. In QED3, therefore, we consider the 
photon propagator of the form: 



kpk v 



D T (k) + ^D L (k) + ep Aj) (k), (10) 



3 



In what follows, / := g implies that / is defined by g. 
4 Note that r anticommutes with 73,75 and commutes with 70,71,72, and that 7] = —75, 7! = — 1. 
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with the non-local gauge [] 



ak 2 



D L (k) 



The vacuum polarization tensor in QED3 is written as 



iV(fc) 



J flV 



k^k v 



n T (k) + ^H L (k) + e^ p ^=Uo(k). 



From the SD equation for the photon propagator, 

[D^k)]- 1 = [D%(k)]- 1 -U^k), 
the general form for Dxih) and Do{k) is written as 

k 2 - U T (k) 



D T (k) 



(k 2 - U T (k)) 2 + (Uo(k) - eVk 2 ) 



(11) 



(12) 



(13) 



(14) 



D {k) 



U (k) - 9Vk 2 



(k 2 - n T (k)) 2 + (Uo(k) - eVk 2 ) 2 



(15) 



Now we define the parity eigenmatrix decomposition by x± = a(l i r )- Then 
we have - B = A e fi + A q t]/> - B e - B q t = (A + fi - B + )x+ + {A-ft - BJ)x~ and 
[A-ft — B]~ l = [A + ft — B + ]~ l x+ + [A-'ft— B-]~ 1 x-- Then the parity even and odd parts 
are obtained as A e = {A + + A_)/2 and A Q = (A + — A-)/2 and similarly for B e and 
B Q : B e = [B + + B_)/2 and B Q = {B + — B_)/2. Thus the SD equation for the fermion 
propagator is decomposed into the simultaneous integral equation for A and B: 



A±(p) 



1 + s 



d 3 q 



1 



p 2 J (2<K)3q 2 A 2 ± (q)+Bi(q) 



±2B ±{q f- 1 ^D (k) 



■P\ 



+A±(q) 



( P -q)t{k 2 )+2 {k - q) ^ P \ l-ak 2 )) 



D T (k)\, (16) 



B± (p) = m± + e 2 J 



d 3 c 



(2^ q 2 A 2 ± (q) + B 2 ± (q) 



B±(q)[2 + i(k 2 )]D T (k) T 2A±(q) q ' {q f } D Q (k) \ , (17) 



\Q — P\ 

where m± := m e db m Q and k := q — p. 

As shown in Appendix A, if we take the following non-local gauge function £(s :- 
k 2 ): 



s 2 D T {s) Jo 



dtD T (Vt)t, 



(18) 



5 This non-local gauge function £ reduces to the usual gauge fixing parameter a of the covariant gauge 
±(dM 2 in the limit D T (k) -» £ and D L {k) -» 
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then we have 



e 



d 3 q 2B ± (q) _ p. (q-p) 



A ± {p) = 1 ± - / \ D (fc) , F ' (19) 

p 2 J (2ity q 2 A z ±(q) + B^q) \q-p\ 

If we do not take into account the CS term Do(k) = 0, the above expression for the 
non-local gauge gives exactly A± = 1 for any B. In the case of the induced CS term, 
i.e., 6 = (see section 4), the term B±(p)Do{k) is of the order B 2 , which is negligible 
in the linearized SD equation which we are going to study in what follows. Even in 
the presence of bare CS term (see section 3), the contribution from the second term 
in the right hand side of eq. ( |i~9| ) is very small compared with 1 in the neighborhood 
of the critical point of the chiral phase transition which we pay special attention to in 
this paper. [] 

If we adopt this non-local gauge, therefore, we obtain A e (p) = 1 and A Q (p) = 0, 
and the even and odd mass functions obey 



B e (p) = m e + e 2 J 



d?q 



(2tt)3 [q 2 + B 2 (q) + B 2 (q)] 2 - AB 2 (q)B 2 {q) 
x j[2 +ak 2 )][q 2 + Bl(q) - B 2 (q)]B e (q)D T (k) 

-AB e (q)B (q) q ' {q ~ P) D (k)\, (20) 



-P\ 



B (p) = m + e 2 J 



d 3 q 



(2vr)3 [q 2 + B 2 (q) + B 2 (q)] 2 - 4B 2 (q)B 2 (q) 
x j[2 +Z(k 2 )}[q 2 ~ B 2 e (q) + B 2 {q))B {q)D T {k) 

+2[q 2 + B 2 e {q) + B 2 (q)] ~ f D (k) j , (21) 

In the non-local gauge given above, B±(p) obey 

B±(p) = m ± + ^J^ dg g2+ ^2 (g) l B ± («)^(P» 1) =F L & ?)] , ( 22 ) 
where we have introduced the ultraviolet (UV) cutoff [] A and defined 



K(p, q) 

and 



/•7T 

/ dd sin 


\ . i{q-p)~ 


Jo 


2 



D T (q-p), (23) 



L(p,q) := T dtisin# q '( q y D (q-p), (24) 



6 This has been confirmed numerically, see p5[ . 

7 It is shown that the solution damps very quickly in the region p > a and hence the UV cutoff can be 



substantially identified with a [jl|, g]. 
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for the angle $ denned by (q — p) 2 = q 2 + p 2 — 2pq cos i?. 

This integral equation is non-linear in B± and difficult to be solved analytically. 
Therefore we adopt the approximation to get the appropriate linear equation. ^| Here 
it should be remarked that, in the linearized equation, if B(p) is a solution, the rescaled 
function nB(p) by a constant k is also a solution when L(p, q) = and m± = 0. Hence 
the magnitude of the linearized integral equation can not be determined in this case. 
Therefore, in the linearized equation, we must introduce an additional condition, the 
normalization condition, which specifies the magnitude of the solution. 

According to [p8| , |29| , we linearize eq. (^) by replacing the denominator q 2 + 
B 2 (q) with q 2 and introduce the infrared (IR) cutoff M in the region of integration: 
M < q < A. Then we adopt the normalization condition M = B(p = M). This 
procedure enforces us to consider the solution only in the region: M < p < A. In this 
approximation, the SD equation reduces to 

B± (p) = m± + £^ £ dq [B± (q)K(p, q) *L(p,q)], (25) 

supplemented with the normalization condition M± = B±{p = M±). 



3 Chiral symmetry breaking 



First of all, we consider the dynamical breaking of the chiral symmetry. In this section, 
we neglect the induced CS term, Tlo(k) = and take into account the one-loop vacuum 
polarization: 



U T (k) = -aVk 2 + 0(B 2 ), a :-- 



e 2 N f 



(26) 



Then D T (k) and D Q (k) read 

j—. k 2 + a\> r k 2 

D T (k) 



{k 2 + aVk 2 ) 2 + 9 2 k 2 



,D (k) 



evk 2 



(k 2 + aVk 2 ) 2 + 9 2 k 2 



(27) 



where k 2 = (p — q) 2 ■= s denotes the squared photon momentum. In this case, eq. (18) 
leads to the non-local gauge function: 



(S + Otyfs) 2 + 6 2 S 
S 2 (s + CHyfs) 



2ay/s + s - Aa9 



a a + y/s 
arctan — — arctan 



(28) 



The function ^ = £(A; 2 ) of k 2 is not singular even at k 2 = and decreases to zero as 
k 2 — ► oo, see Fig.2. The non-local gauge just obtained eq. (p§D has quite complicated 



8 In the linearization of eq. ( J20| ) we should keep the coupling between B e (q) and B (q), which leads to an 
infrared cutoff depending on both B e (0) and B o (0). In the linearization of eq. leading to eq. (p5|), the 
equations for B + and -B_ decouple completely, but it should be noted that M + and M_ are different from 
each other due to the presence of 9 |2q] . 
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form. In order to solve the SD equation analytically, therefore, we expand the non-local 
gauge eq. (p^) around k 2 = and take into account the first and the second terms, 
since it changes rapidly in the infrared (IR) region k/A <C 1. The effect of higher order 
terms will be discussed in the final section. Of course, it is necessary to perform the 
numerical calculation to include the non-local gauge completely, which is a subject of 
the subsequent paper p5| . It is easy to show that £(fc 2 ) of eq. (pg|) has the following 
expansion in the IR region k <C a: 



£(* 2 ) 



1 a 2 - e 2 VW 

3 a 2 + 6 2 a 



+ o 




(29) 



This expansion can be matched with the 1/Nf expansion (when = 0). If we expand 
the non-local gauge into such a series and truncate the series up to 0((k/a) N ), (N > 1), 
the truncated non-local gauge £ guarantees that A(p) = 1 + 0((l/Nf) N+1 ). This 
implies that the WT identity is satisfied at least up to 0((l/Nf) N+1 ) by the vertex 
function T^p, q) = 7^ + 0((l/Nf) N+1 ) and the wavefunction renormalization function 
A{p) = 1 + 0({1/N f) N+1 ). The first term 2/3 corresponds to the gauge choice of 
Carena, Clark and Wagner pif , which is obtained as a consequence of our systematic 
treatment of non-local gauge. 

We substitute the above form eq. (|2^) of the non-local gauge into the SD equation 
eq. ([17|) for B±. In order to solve the equation analytically, we linearize the equation 
based on the second approximation explained in the previous section: 



4a r A 

B±(p) = m± ± / (p + ^— / dqB±(q)K(p, q), 

7T z Nf Jm± 



with 



4o 



ir*l\f Jm± 



a 



dqL{p,q). 



Here K(p, q) is obtained by using eq. ( |29| ) as 



K(p,< 



a 



3a 2 + 6 2 



p + q — \p — q\ 9 q 2 



2pq 



. a: 



+ e 2 



a 



(30) 



(31) 



(32) 



where, corresponding to eq. (|29|) of the non-local gauge-function, we have used the 
following expansion: 

. / 1 \ a 1 a 2 — 9 2 /—t. 



a 2 + Q2 ^ ( a 2 + ^2)2 

Similarly L(p, q) is obtained as 



(33) 



L(p, q) 



a 2 + 6 2 
2a9 



_ p 2 - q 2 p + g 



2pq \p-q\ 



-{p + q-\p 
p 



(a 2 + 6 2 ) 2 

(p + q)(p 2 -pq + q 2 ) - \p - q\ (p 2 + pq + q 2 ) 
3pq 



(34) 



where we have used 

D (k) 



1 



+ 



2a9 



a 2 + 02^2 (a 2 +9 2 ) 2 

Thus we obtain the SD equation 



+ 0(Vk 2 ). 



B±(p)=m±±f(p)+\ dqB ± {q) 

JM± 



p + q — \p — q\ 9 q 2 



1 



2pq 



8a 2 + 9 2 a 



with 



A 



32 



3n 2 N f a 2 + 9 2 ' 



(35) 



(36) 



(37) 



Note that, if 9 = 0, f{p) = and B±(p) obeys the same equation. 

The integral equation eq. ( |36|) can be solved as the boundary value problem of the 
second order differential equation as follows. It can be converted into the differential 
equation 



where 



|( p «) +AB±w . ±rW , 



d ( 2 df(p)\ 

r(p) ■= -r [P — 7— 
dp V dp J 



(38) 



(39) 



Then the solution of the SD integral equation can be obtained from the general solution 
of the differential equation eq. ( |38| ) which satisfies both the IR boundary condition 
(BC): 



B' ± (M±) = ±f'(M±), 



and the UV BC: 



m± = B±(A) + [1 - J(a, 9)] A[B' ± (A) T /'(A)] T /(A), 



where 



9 a 2 
J (a, 9) :=- 



9 2 A 



a 2 + 9 2 a' 



(40) 



(41) 



(42) 



Here it should be remarked that the second term in the non-local gauge function 
eq. ( |29| ) does not change the differential equation eq. ( |38|) and simply changes the UV 
BC. Therefore the second term takes no effect on the general solution of the differential 
equation, but may change the scaling law through the UV BC. Without it, UV BC 
reads m± = #±(A) + A[^(A) T /'(A)] T /(A). 

If r(p) = 0, the general solution of eq. ( |38| ) is given by the linear combination of 
two special solutions: p( _1 + T )/ 2 and p(~ 1 ~ r ^^where r := \/l — 4A. In the presence of 
the inhomogeneous term r(p), the general solution of eq. (|38|) is given by 



B±(p) = M± 



ci 



+ C 2 



±be( P ), 



(43) 
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where 



hip) 



-1 + T 

P 2 



dq q 



-l-T 



1 -1 + T 

dkr(k)k 2 , 



(44) 



Two constants c\ and C2 are determined from the IR BC and the normalization condi- 
tion in what follows. 

In order to write bg(p) explicitly, we must calculate the integration as 



/GO 



3X6 
8a Jo 



dq 



1 



2pq 



In 



p + q 



\P- 



(45) 



Numerical integration shows that the function f(p) is well approximated in the interval 
[0, A] by the function: 

3X9 

m = "8o7 



a A +b{ A ) 



(46) 



with a = 2.5329 and b = 0.7913, see Fig. 3. Especially, in the region p <C A, f(p) 
is well approximated by the linear function with a = 2, b = 0. This can be derived 
by expanding the integrand of f(p) in powers of q/p or p/q after dividing the integral 
region into p > q and p < q, respectively. 
Hence the function r(p) reads 



r{p) 



3X6 
4a 



A 



This leads to 



. , 3X6 
beip) = 



a p 



+ 



3b 



-V 



(47) 



(48) 



2 + AA 6 + AVA/ 

Note that bg (p) is monotonically decreasing in p and bg (p) [ as 6 [ uniformly in p. 
By using f{M) = ^ (-§ + bf ), IR BC eq. © reads 



ci - 



-1 + T 



+ C 2 - 



-1 -r 



3A£ 

" 4a 



X 



+ b- 



X M± 



2(2 + A) 6 + A A 



From the normalization condition, B±(M±) = M±, we obtain 



ci + c 2 = 1 ± 



3X6 

4a 



1 



3 M± 



2 + A 6 + A A 



(49) 



(50) 



Therefore two coefficients c\ , C2 are determined as 



C2 



1 + r 3A6» 



2r 



4ar 



1 - A + r , M± -3 + 2A - 3r 
+ 6- 



-1 + T 3A6< 
4ar 



2r 



2(2 + A) 
a 



A 2(6 + A) 



1 + A + r , M± 3 - 2A - 3r 
+ 



2(2 + A) 



Finally the UV BC eq. © is written as 



[1 + J + t(1 — J)]- 



ci /^M ± \ — 



A 2(6 + A) 

ci / M±\ — 



(51) 



+ [1 + J — r(l — J)] 



A J 



where 



C(X,6) := 1 



A(2-J) fc A(3-2J) 
2(2 + A) 2(6 + A) ' 



(52) 
(53) 
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3.1 # = 



In the absence of 8, B + and B- obey the same SD equation when m± = 0. Therefore 
B e (p) = B±(p) and B a (p) = 0. In this case B±(p) = is a trivial solution of eq. (|17|). 
Here we define M e := M±. Then the UV BC eq. (||) yields 



(54) 



= [1 + J + r(l - J)f-k (%) V + [1 + J ~ r(l - J)]| f ^ " 



A 



when m± = 0. By using the coefficient obtained in eq. (|5l|) : 



A 



ci 



1+T 

2r : 



C2 



-1+T 

2r : 



the ratio M e /A is obtained as 



M P 



A 



(1 + r) [1 + J+(l- J)t] 
(1-r) [l + J-(l- J)r]' 



(55) 



(56) 



The right hand side (RHS) is greater than 1, since < r := v 1 — 4A < 1 and J > 0. 
Therefore this equation can not have the solution M e /A < 1 for real r, i.e., A < 1/4. 
In the region A > A c := 1/4, we define 



uJ 



:= V4A — 1, r = iu, A 



32 



3ir 2 N f ' 



Then we look for the solution of the equation 



1 + J - (1 - J)oj 2 + 2iuj 
1 + J - (1 - J)uj 2 - 2iu 

This equation has the solution < M e /A < 1 of 



X 



Me 

A 



exp 



2ra7r - 2(p 



arctan 



2uj 



1 + J - (1 - J)w 2 



(57) 



(58) 



(59) 



In this region the non-trivial solution is obtained from eq. ( |43| ) and eq. (pl|). 

Thus it is shown that, in the chiral limit m e = 0, the non-trivial chiral-symmetry- 
breaking solution exists only for Nf < N'j as 



B e {p) 



M 3/2 s - m ^ N c /Nf _ Un JL + 5 ) 



sin (5) 



(60) 



with 5 = arctan u = arctan ^jN^/Nf — 1. The scaling law for the fermion dynamical 
mass is given by 



Me 

A 



exp 



2wk 



^ N V N f - 1 



arctan ■ 



2^jNf/N f - 1 



(61) 
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with a critical value for the flavour: 



N 



f 



128 
3^2" 



(62) 



This shows that there exists the finite critical number of flavors Nj above which (Nf > 
Nf) the chiral symmetry restores. 

In the absence of CS term, the value Nj = ||f ~ 4.3 that we have just found is 
4/3 times as large as the value 32 /ir 2 of ANW Q and coincides with that of the 1/Nt 
analysis by Nash pH]. PI This result clearly comes from the first term £ = 2/3 of the 
non-local gauge eq. (|29|). The second term of the non-local gauge does change neither 
the critical value nor the type of the scaling law for the fermion mass, the essential 
singularity at Nf, and merely decreases the critical coefficient, i.e., the magnitude of the 



fermion mass, since M e /A ~ K exp —Inixj J 'Nf/Nf — 1 
the neighborhood of Nf. 



K 



exp 



4/(l + 



'1A 

8 a 



111 



3.2 6^0 



In this subsection we consider the region Nf < Nf = H§. By substituting eq. @ 
into the UV BC eq. (|5^), the equation of state is obtained as 



m± 



[R±(\,9) + iI±(\,0)) 



3+iu 



[R±(X,9)-iI±(X,e)] 



(63) 



where we have defined 

R±(\,9) : = 
I±(X,0) : = 



l + J-co 2 (l-J)± 2G[(1 - A)(l + J) - w 2 (l - J)], 
2w{l±6[2-A(l- J)]}, 



with J = J(a, 6) defined by eq. ([42]) and 



e 



3aA 9 
8(2 + A) a' 



(64) 



(65) 



In deriving this equation eq. (p3|), we have neglected the term with the factor b in 
eq. (ID, since b(M/A) 2 < a(M/A). 

When m-t = 0, we look for the solution for 



From the equation of state, (± obeys 



sin(^ - U - ln C± ) = ±|C ( A,^).K ± 1 / 2 ^ G 3/2 , 
2 2 a 



(66) 



(67) 



3 Nash's non-local gauge corresponds to a special case of ours: Di(k)/DT{k) = constant, independently 



of k. 
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where we have defined 



arctan 



R±{\,oy 



and 



K± = R±(\,6) 2 + I ± (\,9f 



(68) 



(69) 



For the non-trivial solution to exist, the above equation must have the solution for 

o < c± < i. 

We start from the limit 9^0. In this limit, RHS of eq. ( |67| ) vanishes. Hence the 
solution of eq. (|67|) is given by 



t± :-- 



LO 



ln£± = nir (n = 1, 2, 



(70) 



Then, in the region Nf < iVS and = 0, M± obeys the scaling law of the essential 
singularity type: 



M 4 



Aexp 



2nir — 2(f>± 



>0, 



(71) 



as shown in the previous subsection. From the argument 0, ^7]] using the effective 
potential of Cornwall- Jackiw-Tomboulis [34|, it is shown that the solution with n = 1 
corresponds to the nodeless solution and gives the ground state solution and that 
the solution with n > 2 corresponds to the excited solution with nodes n — 1. The 
magnitude of the solution with nodes n is related to the ground state solution as 



M± 




A 



(72) 



This implies the hierarchy of the mass scale: M_£ =1 > M± =2 > M± =3 > .... 

Next we keep the value of Nf below the critical point Nf and increases 6 > 0. To 
find the solution, we rewrite eq. (|6T 



smt± = H±(t±), H±(t) :-- 



as 



±-C(\,6)luK± 



1/2 



\9 ~30± 

— e w 



exp 



3f 

a; 



(73) 



The solution of eq. (|67j ) is given as the intersection point of the sini curve with the 
curve H±(t), see Fig.4. In the absence of the CS term, 8 = 0, eq. fl67p has (countably) 
infinite number of solutions as shown above. However, in the precense of the CS term, 
6/0, the strongly oscillating solutions with many nodes (n 3> 1) disappear and 
eq. (|67| ) has finite number of solutions, no matter how 9 is small, since |fl±(t)| is a 
rapidly and monotonically increasing function in t and H±(t) — > ±oo as t — > oo. 

For small 9 > 0, H+(t) > and H_{t) < for our choice: a = 2.5 and J given by 
eq. (|42]). Paying attention to the region < t < 2ir, we see in this case that t + (resp. 
i_) has a solution at t = ir — 5 + (9) (resp. t = it + slightly smaller (resp. larger) 

than 7r where <5±(0) > and 5±(6 = 0) = 0, see Fig. 4. Hence we obtain the solution 



M± = Ae 



±2 



exp 



2n7r — 2(f>± ' 



(74) 
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Therefore the solution M?. gets larger value than M^ =0 , while M_ gets smaller one than 
M<? =0 . gy substituting the coefficients eq. (|5T|) into eq. (^3|), the non-trivial solution 
is given by 



B±(p) 



M 



3/2 



sin 



sin 5 



± 26- 



sm 



m In 

2 m M± 



+ e 



sine 



±be(p), 



(75) 



where be (p) is given by eq. (|^) and 5 = arctan u and e = arctan . This solution 
satisfies the normalization condition M± 



i-A- 

B±(M±) as should does. 



As discovered in 35], in the presence of the Chern-Simons term, a first order phase 
transition does occur. This phenomenon can be explained in our scheme as follows. 
Consider the situation that 9 is increased gradually with Nf being kept. We find that 
there is a certain value of 9, say 9 c (Nf), such that at 9 = 9 c (Nf) the nodeless solution 
B_{p) (corresponding to n = 1) does disappear suddenly, since H_{t) does not intersect 
with sini any more in the interval < t < tt for sufficiently large 9. At 9 = 9 C , the 
solution B_{p) cease to exist. For 9 > 9 c (Nf), only the nodeless solution B + (p) exists. 
This discontinuity at 9 = 9 c (N f ) persists for B e = {B + +B^)/2 and B Q = (B + -BJ)/2. 
This discontinuity can be observed as the first order phase transition in the direction 
of 9. By using S(p) = [Aj> - B}^ 1 = [A + ^> - B + \~ l x+ + - the chiral 

order parameter {tpip} is written as 



{i>i>)/Nf 



A ^2 



p dp 



7T^ 



B+(p) 



+ 



B-{p) 



p 2 + B + {p) 2 p 2 + B_{p) 



(76) 



Therefore such a first order transition can be observed as the discontinuity of the chiral 
order parameter. 

The critical line is obtained as follows. Substituting eq. (^) back into eq. (p7|), we 
obtain 



X9 



n 



(-1) 



n+l 



2 K l/2 
s in(5 ± (9))^C(\,9) 



-l e ±3 



exp 



3n7r — 30-i 



Up to 0(9), the critical line in the phase diagram (Nf,9) is given by 

3n7r" 



\9 
— ~ exp 

a 



UJ 



(77) 



(78) 



Above this critical line in the phase diagram, the chiral symmetry restores. This shows 
that the CS term has the effect to decrease the critical number of flavors, Nj(9) < 
Nf(9 = 0). However the chiral-symmetry restoring transition is the first order phase 
transition. 

At typical points in the phase diagram, (6,Nf), which is drawn schematically in 
Fig.5, we find (M e = (M+ + M_)/2, M Q = (M+ - M_)/2) 



A 


Nf 




9 = 0;M e = 0,M o 


= 0, M+ = 0, M_ = 0, 


B 


Nf 


<Nf,i 


9 = 0; M e > 0, M 


= 0, M+ = M_ > 0, 


C 


Nf 


<N c f ,t 


9^0;M e >0,M o 


> 0, M+ > M_ > 0, 


D 


Nf 


<N^,i 


9 = 9 c (N f );(first 


order transition) 


E 


Nf 


>Nf,i 


?^0;M e = 0,M o 


= 2M+ > 0, M + = -M- 



> o, 



(79) 
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where the points A, E are in the chiral symmetric phase and B,C,D are in the chiral- 
symmetry-breaking phase. In the above we have considered the route starting at 9 = 
(A —* B — > C — > D) to justify the linearization procedure adopted in this paper, 
since B±(p) = is a trivial solution at 9 = 0. It should be remarked that we can 



adopt another linearization [27] by replacing the denominator q 2 + B 2 (q) with q 2 + M 2 



where M determines the normalization of the solution of the linear SD equation, e.g., 



M = B(0). In [35] and the subsequent paper we give the thorough analysis on the 
choice of linearization procedure which is compatible with the route of approaching 
the critical point, especially, A — > E — > D starting at 9 ^ from the beggining where 
B±(p) = is not a trivial solution. 



4 Parity breaking 



We consider a problem whether or not the parity is dynamically broken in (2+1)- 
dimensional QED, in other words, the non-perturbatively induced CS term may or 
may not break the parity. When 9 = 0, we obtain 



k 2 + aVk 2 
T( ' ~ (^ToVPi^TTM^' 

{k 2 + aVk 2 ) 2 + U (k) 2 

Here it should be remarked that the parity odd vacuum polarization IIo(/c) is at least 
linearly proportional to B a (p), since 

no(fc) = 64 tp/ (0 ( * 2 + "- k) w$h + 0(B ° 2) - <82> 



Therefore the quantity T^j(k) in the denominator of Dj- and Do should be neglected 
in the linearized SD equation. This is sufficient to study the critical phenomena in the 
neighborhood of the critical point. 

From eq. (|2l|), the linearized SD equation for the odd part B Q reads 



^ f ^ q - {q - p) f D (q- P ), (83) 



Nf J (2tt)3 q 2 V{q^pJ 

where the non-local gauge is given by £,(k 2 ) = | — + 0(k 2 ). By mimicking the 
procedure of CCW |HJ], the linearized SD equation can be rewritten as 

32 r A 

B (p) =m Q + J dqK(p, q)B (q), (84) 
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with the kernel 



K(p,q) 



P + Q-\P- 
2pq 

k 2 - 



1 + 



In 



6-7T 

8a 

p + k 



A 1 



a 



k + a 



2pk \p — k\ 



1 + 



k 2 



2qk 



■ m - 



k\ 



(85) 



This integral equation is too difficult to be solved analytically. However it is quite easy 
to see that the SD equation can not have the nodeless (ground state) solution with a 
definite sign, since the kernel K(p, q) is always negative for any p,q £ [0, A] as shown 
in Fig. 6. This implies that the SD equation has only a trivial solution B {p) = 0. 
Therefore the dynamical breakdown of parity does not occur, which is in agreement 
with the previous analyses pi], 22. 23, 24]. 



5 Conclusion and Discussion 

In this paper we have reexamined the dynamical breakdown of chirality and parity in 
(2+l)-dimensional QED. In order to get rid of the ambiguity coming from the vertex 
correction, we have derived the non-local gauge fixing function which guarantees the 
absence of the wavefunction renormalization for the fermion. Using the non-local gauge, 
we have written down explicitly the SD equation for the fermion mass function. 

In the analyses done in this paper, we have taken into account the leading and 
next-to-leading terms of the non-local gauge function £,(k 2 ) in the IR region to solve 
the SD equation analytically. Within this approximation we have found that there 
exists the finite critical number of flavors Nj = 128/37T 2 and that the scaling form of 
the essential singularity type is obtained in the absence of the CS term. The critical 



value N'j obtained in this paper reproduces Nash's result [3C] which is slightly larger 
than the previous one 32/tt 2 in the absence of the CS term 6 = 0. 



In the presence of the CS term, it is shown [35] that the phase transition in the 
direction of 9 shows the discontinuity at some critical value, 6 c (Nf), which depends 
on Nf. The bare CS term decreases the critical flavour number. The order of chiral 
transition as well as the critical line obtained in fBf| is completely different from that of 



Hong and Park [33]. In Hong and Park [33] the first order term of 6 was neglected in the 
SD equation for the fermion mass function which corresponds to eq. ( ^o|) and hence this 
approximation inevitably forces the parity-odd mass part to decouple completely from 
the SD equation for the parity-even mass. This careless treatment leads to the totally 
different conclusion: 1) the critical number of flavors as a function of 0: Nf(0) = 

Nf l+ (g/ a y2 , which implies the critical line £ = yjNJ/Nf — 1 where Nj = Nj{9 = 

0). 2) the continuous (infinite order) phase transition characterized by the essential 

singularity on the whole critical line: ^ ~ exp — j==== = , as in the absence of 

CS term. However this treatment is not correct as was noticed in the footnote of the 
linear approximation in section 2, since the CS term dominates at large momentum 
and gives the essential contribution to the UVBC which determines the character of 
the phase transition. This fact was confirmed by the direct numerical calculation of 
the non- linear SD equation [35| without relying on the specific linearization procedure. 
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In the same scheme we have shown that the dynamical breakdown of parity does 



not occur, which agrees with the previous analyses [21, ^2], 23, 24] and the general 



argument [25]. An interpretation of absence of dynamical parity breaking is as follows. 
The gauge interaction acts as an attractive force to cause the (chiral) condensation 
between fermion and antifermion, (tpip) ^ 0. As can be seen from the analysis in 
the previous section, however, the induced CS term generates the same effect as the 
repulsive force and this gets superior to the original attractive force to result in the net 
repulsive force (minus sign of the kernel) which prevents the condensation, {ipTip) = 0. 
Therefore it is necessary to incorporate other attractive interactions so as to induce 
the parity breakdown (ipTip) ^ 0, e.g., the four-fermion interaction as done in CCW 



m 



We must comment on the limitation of our approach. The explicit form of the 
non-local gauge we have found in this paper can not be applied to the vertex which is 
not of the form T^(p, q) = 7„.F(£>, q) where F(p, q) is the functional of A and has the 
limit: F(p,q) — > 1 as A(p) — > 1. For the given vertex which has the different tensor 
structure |i~l|| , therefore, we must rederive the corresponding non-local gauge, if this is 
possible. In such a case, we can not say anything about the result at this stage. 

In this paper we have truncated the 1/Nf series up to some finite order and have 
taken into account the first two terms. To improve our analyses, we should mention 
the effect of the third term in the non-local gauge: 

2 . 2 la 2 -6 2 V& , (a 2 — 59 2 )(3a 2 + 6 2 ) k 2 ( ( k\ 3 ^ 



^ ) = 3-3aTT^T + lHa 2 + e 2 ) 2 a 2+ °{{a) J" (§6) 

The third term changes the equivalent differential equation and hence the solution of 
the SD equation may change. It may be expected that the qualitative feature of the 
dynamical symmetry breaking will be unchanged. If the 1/Nf series can be summed 
as pointed out in |ll|| , however, the wavefunction renormalization effect can alter the 
behavior of the dynamically important low momentum region, which is recognized also 



in [1C[. Moreover we have neglected all the terms of the order 0(B±) to obtain the 
closed SD equation for B e and B Q . The resulting equation is a linear integral equation 
where the IR cutoff must be introduced to define the equation in the IR region. This 
procedure may change the low momentum behavior of the solution, in contrast to the 
UV region. Anyway it is quite necessary to perform the numerical calculation based 
on the non-local gauge without the expansion in order to confirm our results. 

Finally we mention another problem. Though the inclusion of the bare CS term 
breaks the parity explicitly, it is an interesting question whether the induced CS term 
recovers the parity invariance which is explicitly broken at the tree level, 6^0. In 
order to answer this question, we must include both the bare and the induced CS 
terms, in which case it is difficult to solve the corresponding SD equation analytically. 
This problem will be discussed in a forthcoming paper. 
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Research from the Ministry of Education, Science and Culture (No. 05640339), and 
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comments on the effect of the Chern-Simons term. 
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A Derivation of the non-local gauge 

In QED, we introduce the non-local gauge fixing function £(fc 2 ) in the photon propa- 
gator: 



where 



d(k 2 ) := k 2 D T (k), V (k 2 ) := 1 - £(k 2 ). 



(87) 



(88) 



In the euclidean QED/) in D-dimensional space, the SD equation for the fermion wave 
function renormalization A is given by 



P 2 A( P 2 )-p 2 
d D q 



f d g 
J (2tt) d 



A(q 2 



(2ir) D q 2 A 2 (q 2 ) + B 2 (q 2 ) 



xd(k 2 ) 



(D-2)^ + 



p-q (p-q p 2 q 2 - (p ■ q) 



- 2- 



k 2 fc 4 
and the SD equation for the fermion mass function B reads 



V (k 2 ) 



B(p 2 



' f dD( l 
J (2tt) d 



B{q 2 ) 



d(k 2 



■[D-r](k 2 



(2ir) D q 2 A 2 {q 2 ) + B 2 {q 2 ) k 2 
Separating the angle ■& defined by 

k 2 := (q - p) 2 = x + y - 2^/xy cosi?, x:=p 2 , y:=q 2 , 



we find 



p 2 A(p 2 ) - p 2 
rA 2 

C D e 2 / dy 
Jo 

xd(k 2 ) 



2)/2 A{y) 



Jo 



d$sm D - 2 $ 



yA 2 (y) + B 2 (y) 
~xy cos "d{D — 2 + rj(k 2 )) xy — {*Jxy cost?) 2 



k 2 



k 4 



r](k 2 



and 



B(p 2 ) = C D e 2 



where 



C D :-- 



A2 y^-^^Bjy) d(k 2 ) 
y yA 2 (y)+B 2 (y) k 2 

1 



[D-n(k 2 )}, 



2 D 7r (D+l)/2 T (D=iy 
We perform the angular integration by parts as 



rn 



d"&sm D 2 -d cos -d yjxy f (z) 



xy 



D - 1 



sm D - 1 ^f{z) 



t 2xy 
~ D - 1 



(89) 



(90) 



(91) 



(92) 

(93) 
(94) 



cMsin D tif'(z), 
(95) 
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where 

z :=k 2 = {q-pf. (96) 

Then we find 

2^ f 2^ 2 2^ J V {D - 2),2 A{y) (2xy) r j a ■ D a 

pA{p) - p = - eCD L dy Wrf(#/o Mm d 



1 d f(D -2 + r](z))d(z)\ d(z)rj(z) 



D — 1 dz V z J z 



-.2 



(97) 



This is further rewritten as 



2 A( 2, 2 C D e 2 /"A 2 y( D -^A(y)(2xy) r . D 

P ^ P ^~ P = -D^iJa dy y A 2 (y) + B 2 (y) L * 



X ^=T {(z D - 2 d(z)r,(z))' — (D — 2)z D -*[d(z) - zd'(z)}} , (98) 

where the prime denotes the differentiation with respect to z. 

The requirement A(p 2 ) = 1 is achieved by taking rj(z) such that {z D ~ 2 d{z)r){z))' — 
(D — 2)z D ~^[d{z) — zd'(z)] = 0. This is simply solved as follows. 

V{*) = oSm , r dt[d{t)-td'{t)\t D -\ (99) 



z D - 2 d{z) Jo 

where we have assumed that [z D ~ 2 d(z)r](z)]\ z= o = so as to eliminate the l/z D ~ 2 
singularity in rj(z) and £(z) = 1 — rj(z). This should be checked after having obtained 
the function rj(z). 

The gauge-fixing function = 1 — rj(z) is simplified as 

( D- 1)(D- 2) f* MJUuD _ 3 
z D ~ 2 d{z) 

where we have assumed [z D ~ 2 d(z)]\ z= o = 0. 

For D = 3, eq. fllOCl ) reduces to eq. @), since k*D T {k 2 ), k A D T {k 2 )£{k 2 ) -> as 
k 2 ^0. 



£( Z ) = D-l- y g _" > \ dtd{t)t D ~\ (100) 



B Integration formulae 

Define the angle i? from k 2 := (q — p) 2 = p 2 + q 2 — 2pqcos , d, where k = \^k^,p 
\fjfi", q = \/?- Then we obtain the following results. 

. ,1 p + q-\p-q\ 2 



ok pq max(p,q) 

cos i? (p + q)(p 2 - pq + q 2 ) — \p— q\(p 2 + pq + q 2 ) 



dfismfi- 



o 



k 3p 2 q 2 
2 min(p, q) 



,2 „2\ 



3 max(p 2 , 
ausmu-rr = — in- 



o 



A; 2 |p — g| 



' • qCOST? 1 P 2 + ^ 2 P + q . . 

dvsmv—z- = h r. 9 o m l r- 

fc z 2jrg z |p — <jr| 
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Figure Captions 

Fig.l: Schwinger-Dyson equation for the fermion propagator. 

Fig. 2: Non-local gauge function £(s) for s := k 2 (when 9 = 0). The lines from above to 
below correspond to a/A = 1, 10 _1 , 1CT 2 , 10~ 4 . 

Fig. 3: Plot of the function f(p) in eq. (fj5|), apart from the factor of the integral . 

Fig. 4: Schematic plot of the function H±(t) eq. ([73]) and sini. The solution M± is 
obtained from the intersection point of the curve H±(t) with sini. 

Fig. 5: Phase diagram (schematic). 

Fig. 6: Integral kernel K(p, q) as a function of p and q where A = a. 
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